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ABSTRACT 

Based on cosmological rates, it is probable that at least once in the last Gy 
the Earth has been irradiated by a gamma-ray burst in our Galaxy from within 2 
kpc. Using a two-dimensional atmospheric model we have computed the effects 
upon the Earth’s atmosphere of one such burst. A ten second burst delivering 
100 kJ/m 2 to the Earth results in globally averaged ozone depletion of 35%, with 
depletion reaching 55% at some latitudes. Significant global depletion persists 
for over 5 years after the burst. This depletion would have dramatic implications 
for life since a 50% decrease in ozone column density results in approximately 
three times the normal UVB flux. Widespread extinctions are likely, based on 
extrapolation from UVB sensitivity of modern organisms. 

Subject headings: astrobiology - gamma rays: bursts 


1. Introduction 

Camma-ray bursts (GRB) within our galaxy have been suggested as a possible threat 
to life on Earth (Thorsett 1995; Scalo & Wheeler 2002; Dar &; De Rujula 2002; Melott et al. 
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2004). Some effects similar to that due to a nearby supernova (SN) are expected, especially 
depletion of stratospheric ozone due to ionization caused by incident gamma radiation. GRB 
are rarer than SNe, but their greater energy output results in a larger region of influence 
and hence they may pose a greater threat. Results presented in Melott et al. (2004) indicate 
that a GRB impact is likely within the last Gy. Such a GRB is expected to have dramatic 
effects on stratospheric ozone, leading to detrimental effects on life through increases in solar 
ultraviolet (UV) radiation which is strongly absorbed by ozone. 

In order to gain more detailed and accurate insight into these expected effects we have 
performed computations using the Goddard Space Flight Center (GSFC) two-dimensional 
atmospheric model. This model has been used previously to investigate the atmospheric 
effects of SN (Gehrels et al. 2003) and we build on that work. 


2. Methods 

We take as “typical” a GRB with power 5x 10 44 W (isotropic equivalent) and duration 10 
s, whose gamma-ray spectrum is described by the Band spectrum (Band et al. 1993). These 
assumptions are drawn from observations and are not dependent upon beaming angle. We 
have determined the depletion of ozone for such a GRB beamed at the Earth from a distance 
of 2 kpc, delivering to the Earth a total fluence of 100 kJ/m 2 . This distance corresponds to 
that of a probable nearest “typical” GRB in the last Gy, based on conservative assumptions 
(Melott et al. 2004). 

The prompt effect of this burst at the Earth’s surface is a “flash” of UVB radiation 
with power ^ 20 W/m 2 (Smith et al. 2003). This is about seven times the intensity at the 
Earth’s surface on a bright, sunny day, but is brief and so is not likely to have a major effect 
on life. Longer term effects include ozone depletion and the resulting increase in solar UVB 
flux, which we begin to explore here. 

We have not included the effects of any ultra high energy (> 10 18 eV) cosmic rays from 
a GRB (Dermer & Atoyan 2004), due to uncertainty as to whether and at what energies GRB 
may produce such particles. Lower energy cosmic rays which may be produced are negligible 
at our assumed distances due to deflection by the galactic magnetic field. It is estimated 
(Dermer 2004) that inclusion of ultra high energy cosmic rays would add about 30% to the 
ozone-destroying energy budget. We would also expect production of some radioisotopes 
from cosmic ray spallation off of atmospheric constituents, as well as a terrestrial burst of 


muons. 
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2.1. Atmospheric Model 

The GSFC 2D model is described in Douglass, Jackman, & Stolarski (1989); Jackman 
et al. (1990); and Considine, Douglass, & Jackman (1994). The model’s two dimensions 
are latitude and altitude (ranging up to about 116 km). The latitude range is divided 
into 18 equal bands and extends from pole to pole. The altitude range includes 58 evenly 
spaced logarithmic pressure levels (approximately 2 km spacing). A lookup table is used 
for computation of . photolytic source term, used in calculations of photodissociation rates 
of atmospheric constituents by sunlight (Jackman et al. 1996). Winds and small scale mix- 
ing are included as described in Fleming et al. (1999). For this study, we have removed 
anthropogenic compounds such as CFCs. 

We have employed two versions of the atmospheric model. One is intended for long 
term runs (many years) and includes all transport mechanisms (e.g., winds and diffusion). 
This version has a time step of one day and computes daily averaged constituent values. 
The second version is used for short term runs (a few days) and calculates constituent values 
throughout the day and night, but does not include transport. Previously, this version has 
been used with a time step of 225 seconds (Jackman et al. 2001). In the current study we 
have used a time step of one second in order to allow for inputting our GRB gamma radiation 
over several time steps. 

Gamma-rays are introduced in the model in a manner similar to that described in 
Gehrels et al. (2003). In that study, gamma-rays were included using the spectrum of SN 
1987A. In the present study, the gamma-ray differential photon count spectrum used is that 
of Band et al. (1993), which consists of two smoothly connected power laws. We use the 
following values for the break energy and power law indices, respectively: E 0 = 187.5 keV, 
a ~ —0.8, /? =■ —2.3 (Preece et al. 2000). The total received energy is scaled to our desired 
value (in this study, corresponding to a fluence of 100 kJ/m 2 ) The total photon flux in each 
of 66 evenly spaced logarithmic energy bins, ranging 0.001 ^ E ^ 10 MeV, is obtained by 
integrating the Band spectrum for each bin. 


2.2. Simulations 

All simulation runs used for analysis were begun with initial conditions obtained from 
a long-term (roughly 40 years) run intended to bring the model to equilibrium. Constituent 
values from this run are read in by the Is time step version of the model which runs for 7 days 
(beginning at noon), either with or without input of gamma-rays. Runs including gamma 
radiation input the burst as a step function at noon on day 4, with duration 10 s. In the 
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current study, the burst is input in late March (near the spring equinox) with incident angle 
0 degrees (equatorial). (Forthcoming studies will investigate the effects of varying intensity, 
incidence angle, and times of year at which the burst occurs.) Day 4 is the middle of the 
short term run, which allows for “warm up” in the model and insures that relevant chemistry 
is accurately computed over many time steps after the burst. Constituent values from this 
type of run are then read in by the 1 day time step version which is run for 20 years in order 
to investigate long term effects and determine how long the atmosphere takes to return to 
equilibrium, pre-burst conditions. 

Analysis is then performed by comparing such a combined base-short-long run without 
gamma-ray input to such a run with the burst included. Ozone depletion and other effects 
are computed by comparing these two runs. 


3. Results 

Stratospheric ozone is lost through several catalytic reactions involving oxygen-, nitrogen- 
, hydrogen-, chlorine-, and bromine-containing gases. The constituents in the stratosphere 
are generally grouped into “families” such as O x (O3, O, 0( 1 D)), NO y (N, NO, NO2, NO3, 
N 2 Os, HNO3, HO2NO2, CIONO2, Br0N02), HO x (H, OH, HO2), Cly (chlorine-containing 
inorganic molecules), and Bry (bromine-containing inorganic molecules), which allow for ef- 
ficient computation of the chemistry and transport effects. We assumed that there were 
no human-made sources for any of these families, as such large-scale anthropogenic per- 
turbations have been present only over a relatively brief period of time (« decades). Our 
model simulations are thus focused on. a period of geologic time wherein the ozone layer was 
substantial, but the human influence was negligible. 

In the case of a large input of gamma rays to the atmosphere NO y compounds (most 
importantly NO and NO2) are created through dissociation of N 2 which then reacts with 0 2 
to generate NO. Subsequent reactions create NO2 and other compounds. Together, these 
react catalytically to deplete O3 through the cycle 


NO + O3 —y NO2 H - O2 

(i) 

NO2 + 0 — y NO O2 

(2) 

net : O3 + 0 O2 + O2 

( 3 ) 


Note that NO is not consumed in this cycle and the net result is the destruction of O3 
and production of O2. Other reactions can complicate this cycle, such as destruction of NO 
by reaction with N; production of O3 through reactions of NO with HO2; and interference of 
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NO y with other families (chlorine-, bromine-, and hydrogen-containing constituents) which 
reduces the ozone depletion from these families. Some uncertainties in the atmospheric 
model’s treatment of this cycle are discussed in Section 4. 

The primary results of our simulations are increases in NO y and decreases in O3. Ozone 
column densities can then be used to calculate the resulting UVB flux at the Earth’s surface 
and from this information biological effects can be estimated. UVB is particularly damaging 
to organisms because DNA absorbs strongly in this wavelength region and is damaged. 

Results of our modeling are shown in Figures 1 and 2. We have modeled the effects of 
100 kJ/m 2 total received gamma-ray fluence, input as described in Sect. 2.2. This corre- 
sponds to our “typical” GRB located at about 2 kpc. Figure 1 shows the percent difference 
(between a run with gamma-ray input and one without) in vertical column density of O3. 
The burst is input at time 0 and immediate depletion of ozone is evident. Due to their 
qualitative similarity, we have chosen not to plot here changes in NO y . Maximum increase 
in NOy is largely coincident with maximum decrease in O3. The maximum increase in NO y 
is approximately 30-fold for this case. 

Figure 2 shows the vertical column density of O3 at each latitude over time. Again, 
the burst is input at time 0. Included are scales in both Dobson units (the usual unit of 
ozone column density) and 10 18 cm -2 . A Dobson unit describes the thickness of a column 
of ozone at standard temperature and pressure and is defined as 1 DU = 0.01 mm thickness 
(or, 1 DU = 2.69 x 10 18 cm“ 2 ). 

Several features in these plots are worth noting. First, as is seen in Fig. 1 and 2, 
depletion of ozone is initially greatest at the equator (where the incident flux is highest), 
becoming greatest toward the poles within a year or so. Larger ozone depletions at the 
poles are primarily due to the long lifetime of the enhanced NO y in the polar stratosphere. 
Fig. 1 gives a somewhat exaggerated impression of the effect of depletion at the poles, 
since ozone is initially high there. The enhanced NO y , including HNO3, will lead to an 
enhancement of nitric acid trihydrate (NAT) polar stratospheric clouds (PSCs). These NAT 
PSCs facilitate heterogeneous reactions that result in greater ozone depletion by halogen 
(chlorine and bromine) constituents. This is especially true in the south polar region where 
a stronger polar vortex with colder stratospheric temperatures is in place during the winter. 
This contributes to ah asymmetry which would be peculiar to the present-day configuration 
of continents. An additional effect which contributes to the polar asymmetry is the time of 
year at which the burst occurs, since ozone concentrations at the poles exhibit large seasonal 
variations. 

Note, however, that around 5-6 months after the burst there is a short term production of 
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ozone toward the south pole. Production occurs at the end of south polar night when a lack 
of photolysis has caused accumulation of N0 y constituents which are suddenly photolysized 
as the sun rises, producing O which may then react with O 2 to form O 3 . This production is 
fairly short lived as the O is consumed by various reactions and remaining NO y compounds 
continue to destroy O 3 . 

Overall, globally averaged ozone depletion reaches about 35% (at the start of the long- 
term run) and a maximum depletion of about 55% occurs first at the equator immediately 
after the burst, and then again about 15 months after the burst, in the southern hemisphere. 
Significant global depletion (10% or more) lasts for over 5 years after the burst. 

Figure 3 shows DNA damage estimated by convolving the daily average UVB flux with 
a biological weighting function (Setlow 1974; Smith et al. 1980). We have normalized the 
plot by dividing the damage by the annual global average damage in the absence of a GRB. 
Greater DNA damage probability is evident at low latitudes. This is due to combination of 
the O 3 depletion effects with the sun incidence angle, length of day, etc. One might think that 
this damage would be countered by a greater evolved UVB-resistance in organisms at low 
latitudes. However, at least for modern phytoplankton, there is no evidence that temperate 
zone phytoplankton are any more UVB-resistant than Antarctic plankton (Prezelin 2004). 
Thus, at least for this burst geometry, one might predict that greater ecological damage and 
extinction would be likely near the equator. Of course, early in the Earth’s history O 2 , and 
hence 0 3 , levels were lower than today (Thomas et al. (2004) and references therein) and 
life may have been more UVB-resistant, lessening this impact. 


4. Uncertainties 

Our ionization profiles are computed- using simple energy-dependent attenuation co- 
efficients, instead of a full radiative transfer calculation. This technique is implemented 
following Gehrels et al. (2003) with the primary modification being the functional form of 
the spectrum. That study found good agreement between their calculations and a full ra- 
diative transfer model. Similarly, we find that ionization due to the gamma-ray input peaks 
around 30 km elevation, which is in agreement with Gehrels et al. and Smith et al. (2003). 

A discussion of some uncertainties in the production of NO y compounds is presented in 
Melott et al. (2004). In particular, reactions involving exited state nitrogen atoms, N( 2 D), are 
not included in the atmospheric model. Comparison with simplified, off-line computations 
using various ratios of N( 2 D)/N( 4 S) indicate that for our present case the effect on NO 
production of including excited state N atoms is small. Though results depend on the 
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assumed temperature at which the reactions occur, for a reasonable range of temperatures 
(230-270 K), there is little variation of NO production with increased concentrations of 
excited state N atoms, and there is generally less than a factor of two difference between the 
off-line computations and model results. 

The GSFC two-dimensional atmospheric model is empirically based and its dynamics 
are not coupled to the significantly changing constituent levels and accompanying heating. 
This fact introduces some uncertainty in the transport of constituents. 

5. Discussion 

A significant result of our modeling is that even for a short duration input of radia- 
tion, atmospheric effects are large and long-lived. Expectations based on supernova studies 
indicated that a short duration input might not have such large effects (Gehrels et al. 2003). 

Appearance of features such as the localized production of O 3 highlights the need for 
detailed modeling of the effects of a GRB on the Earth’s atmosphere, as complexities and 
feedbacks in the chemical processes can be important and simple scaling may not prove 
reliable. 

Melott et al. (2004) summarizes studies of UVB sensitivity of various organisms (see 
also Cockell (1999)). About 90% of UVB is presently absorbed by atmospheric ozone. Due 
to the sensitivity of DNA to this radiation, increases of only 10-30% can have lethal effects 
on many organisms, especially phytoplankton, the base of the food chain. Ozone depletions 
in the range of 50%, as seen here, lead to roughly three times more UVB at the surface, 
which is clearly a possible candidate for causing mass extinctions. 
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Fig. 1. — Percent change in column density of ozone (comparing runs with and without 
burst) . 
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. — Column density of 0 3 with scales for both Dobson units (left) and 10 18 cm 
(Plotted from one year before the burst to four years after.) 
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Fig. 3. — Relative DNA damage (dimensionless), normalized by the annual global average 
damage in the absence of a GRB. (Plotted from one year before the burst to four years 
after.) 
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Brief, Popular Summary of the Paper: 

This paper concerns the possibility that a very large gamma-ray burst (GRB) in our 
Galaxy influenced the Earth’s atmosphere in the past billion years. The GSFC two-dimensional 
(latitude vs. altitude) model was used in this study to investigate this scenario. 

A very large ten second burst of gamma and X-rays was simulated to impact the Earth’s 
atmosphere. This radiation produced huge amounts nitrogen-containing compounds, over an 
order of magnitude larger than background amounts. These nitrogen-containing compounds 
(primarily nitric oxide and nitrogen dioxide) led to huge amounts of ozone destruction in the 
stratosphere. Globally averaged total ozone depletion was computed to be a maximum of about 
35%. Since the nitrogen-containing compounds can have a fairly long lifetime, significant 
global ozone depletion (>10%) was computed to last for over 5 years past the GRB. 

An associated enhancement in ultraviolet (UV) radiation was computed to be one result 
of the reduced ozone. These enhanced levels of UV can harm various organisms possibly 
leading to extinction of some species. This study thus suggests that a very large GRB in the 
Earth’s past may have led to a mass extinction event. 


